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ABSTRACT 

Using a set of deep imaging data obtained by the Advanced Camera for Surveys (ACS) on the 
Hubble Space Telescope (HST) shortly after its deployment, Yan, Windhorst & Cohen (2003) found 
a large number of F775W-band dropouts (i-dropouts), which are consistent with being galaxies at 
z « 6. The surface density of i-dropouts thus derived, however, is an order of magnitude higher 
than those subsequent studies found in other deep ACS fields, including the Hubble Ultra-Deep Field 
(HUDF). Here we revisit this problem, using both the existing data and the new data obtained for 
this purpose. We confirm that the large overdensity of i-dropouts does exist in this field, and that 
their optical-to-IR colors are similar to those in the HUDF. However, we have discovered that the 
i-dropout overdensity is accompanied with an even larger excess of faint field objects in this region and 
its vicinity. This large excess of field objects is most likely caused by the fact that we have resolved 
the faint diffuse light extending from an interacting galaxy pair in the Virgo Cluster, M60/NGC4647, 
which lies several arcminutes away from the region where the excess is found. The integrated light 
from the excess is a few percent of the luminosity of the galaxy pair. We argue that this excess 
population is most likely tidal "debris" and/or halo stars related to the galaxy pair rather than to 
the Virgo Cluster in general. The z-dropouts in this field are within the magnitude range where this 
excess of field objects occurs. The spatial distribution of the i-dropouts seems to follow the same 
gradient as the excess field population. This excess population is also red in color, and the red wing 
of its color distribution continuously extends to the regime where the i-dropouts reside. While we 
still cannot completely rule out the possibility that the overdensity of i-dropouts might be a genuine 
large-scale structure of galaxies at z « 6, we prefer the interpretation that most of them are part of 
the excess stellar population related to M60/NGC4647. Future spectroscopic work will be needed to 
unambiguously identify the nature of this i-dropout overdensity. 

Subject headings: galaxies: high-redshift — galaxies: halos — galaxies: interactions — stars: AGB 
and post- AGB 



1. INTRODUCTION 

The so-called "Lyman-break" technique (Steidel & 
Hamilton 1992) is now widely used to select galaxies at 
high redshifts. This method relies on multi-band imaging 
to identify the characteristic discontinuity — the Lyman- 
break — in the spectral energy distributions (SED) of 
high-redshift galaxies, which is largely caused by the 
Lyman-limit and Lya absorptions of intervening neu- 
tral hydrogen along the sightlines to such galaxies (e.g., 
Madau 1995). The technique requires imaging in at least 
two passbands, one to the blue side of the break and the 
other to the red side. Lyman-break makes high-redshift 
galaxies much fainter in the blue band than in the red 
one, or in other words, it makes them seem to "drop-out" 
from the blue band. For this reason, this method is also 
known as the "drop-out" selection, and the candidates 
found in this way are usually referred to as "dropouts" . 
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As its first application in the z « 6 regime, Yan, Wind- 
horst & Cohen (2003; hereafter YWC03) used this tech- 
nique in a deep field observed by the Advanced Camera 
for Surveys (ACS; Ford et al. 2003) in the default pure 
parallel mode (Sparks et al. 2001) soon after its deploy- 
ment on-board the Hubble Space Telescope {HSTj. In 
this redshift range, Lyman-break occurs at around 8512A 
in observer's frame (restframe 1216A), which is well tar- 
geted by the F775W and F850LP filters through which 
this deep parallel field was imaged. YWC03 found a large 
number of F775W-band dropouts, and argued that the 
vast majority of them were very likely galaxies at z « 6. 

The number density of z « 6 galaxies inferred from 
YWC03, however, seem to be much higher than those in 
other deep ACS fields (e.g., Bouwens et al. 2003; Stan- 
way, Bunker & McMahon 2003; Dickinson et al. 2004; 
Bouwens et al. 2004; Bunker et al. 2004; Yan & Wind- 
horst 2004, hereafter YW04). While the cosmic variance 
could result in significant fluctuation in the number den- 
sity (e.g., Somerville et al. 2004; Bouwens et al. 2004), 
it cannot account for the large difference between the re- 
sult of YWC03 and those of others. As the true number 
density of galaxies at z « 6 is closely related to a series 
of important cosmological questions (e.g., the formation 
of early galaxies; the reionization history of the universe; 
etc.), it is prudent to closely examine the cause of the 
observed difference. 

In this paper, we present new results that shed light 
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to the nature of this ovcrdcnsity. We first describe the 
existing and new HST observations in §2 and the data re- 
duction in §3. The overdensity of F775W-band dropouts 
in YWC03 is scrutinized in §4, and a new interpreta- 
tion is given in §5. We conclude with a summary in 
§6. For simpHcity, we refer to the HST F475W, F775W, 
F850LP, FllOW, and F160W passbands as 5475, ^775, 
•^850 7 •/no, and -ffieo, respectively. The F775W-band 
dropouts are then referred to as z-dropouts. The mag- 
nitudes are in AB system. All coordinates quoted are 
in J2000. The following cosmological parameters from 
the first-year Wilkinson Microwave Anisotropy Probe 
(WMAP) results in Spergal ct al. (2003) arc adopted: 
lnM,^A,Ho) = (0.27,0.73,71). Using the three-year 
WMAP results (Spergal et al. 2007) would not change 
any of our results. 

2. OBSERVATIONS 

In addition to the ACS parallel field discussed in 
YWC03 (hereafter referred to as "Pari"), we have ac- 
quired two sets of new observations to help tackle the 
question at hand. To follow-up a subset of the i-dropouts 
found in YWC03, two regions in Pari were imaged by 
the Near Infrared Camera and Multi-Object Spectrome- 
ter (NICMOS) Camera 3 (NIC3). During the course of 
these NIC3 observations, the ACS instrument was work- 
ing in parallel, which resulted in a new ACS parallel field 
(hereafter referred to as "Par2") that is several arcmin- 
utcs away from the original one. The layout of these two 
ACS fields are show in Fig. 1. Note that they are in the 
general direction of the Virgo Cluster. We describe all 
these observations below. 

2.1. ACS parallel observations of Pari 

The center of Pari is at RA = 12'*43™32% Dec = 
11°40'32". The ACS data were taken during HST Cy- 
cle 11 soon after the instrument was installed in March, 
2002. In fact, the observations spanned from April 28 to 
June 19, 2002, and were acquired during the execution of 
the Guest Observer (GO) program ID 9043 (PI. Tonry), 
when it was using the Wide-Field Planetary Camera 2 
(WFPC2) as the primary instrument to image NGC4647, 
which is a member of the Virgo Cluster. Only 1775 and 
Z850 filters were used for these ACS parallels. In total, 
15 images were taken in «775 and 27 images were taken 
in ^850- The total exposure time in these two bands is 
2.65 hours and 4.28 hours, respectively. 

2.2. NIC3 observatAons in Pari 

Our NIC3 observations of Pari were carried out as a 
GO program (PID 9780) in HST Cycle-12. Two regions 
in Pari, designated as fields "Subl" and "Sub2", were 
observed in Jno and Hiqq. Due to an unfortunate sys- 
tem glitch that caused NICMOS enter a short safing- 
mode period in early August 2003, our program missed 
its originally planned visit windows, and its execution 
had to stretch from January to December 2004 through 
rescheduled visits. 

For each field, the observations were done at 15 dither- 
ing positions in Jno and 25 positions in iJieo- The ob- 
servation mode was MULTIACCUM with SPARS64 se- 
quence. The exposure time at each dithering position is 
512 seconds, resulting in a total integration time of 2.13 
hours in Jno and 3.56 hours in iJieo, respectively. 



2.3. ACS parallel observations of Par2 

While our NIC3 program was not executed in August 

2003 as originally scheduled, part of its ACS parallel ob- 
servations were still automatically carried out on August 
5 and 6, 2003 as planned in default. This resulted in 
Par2, which is about 8' to the north of Pari. The field 
center is RA = 12'*43'"30^ Dec = ll°49'21^ In total, 
5 images in 2:850, 12 images in 1775 and 6 images in 5475 
were taken in this field, giving a net exposure time of 
1.00, 1.80 and 0.89 hours in these three bands, respec- 
tively. 

3. DATA REDUCTION AND PHOTOMETRY 

In this section, we discuss the data reduction and pho- 
tometry of all relevant HST images. 

3.1. AGS data reduction 

When the YWC03 paper was written (September 
2002), the HST "On-the-Fly Reprocessing" (OTFR) cal- 
ibration pipeline could only process ACS data up to the 
step of geometric distortion correction. The Multidriz- 
zle package (Kockemoer et al. 2002), now a standard 
tool for mosaicing ACS images using the drizzle algo- 
rithm (Fruchter & Hook 2002), was then not yet well 
known and tested. Therefore, the ACS mosaics of Pari 
used in YWC03 were created by stacking the flat-fielded, 
distortion-corrected images utilizing the well-tested IM- 
GOMBINE task of IRAF that had been used by the com- 
munity for decades. 

For this current study, the ACS data of Pari were re- 
trieved again from the HST data archive. In so doing, 
the requested data were automatically reprocessed by the 
OTFR pipeline using the best reference files (bias, dark 
current, flat field, etc.) currently available. New mosaics 
were then made from these reprocessed images by using 
the now-fully-matured Multidrizzle routine as available 
in STSDAS ' under PyRAF. The reduction of the Par2 
ACS data followed the same procedure. 

Some ACS/WFC OTFR-processcd images are known 
to show a quadrant-to-quadrant "jump" in background 
because of the varing residual bias levels in the four 
ACS/WFC readout amplifiers (Pavlovsky et al. 2005). 
While this bias offset is minimal in the data of Par2 
(~ 0.1e~ level), it is obvious in most of the Pari im- 
ages (~ 10e~ level). Although such "jumps" do not 
have significant impact on the photometry of compact 
sources in general, they do affect the sources that are 
close to quadrant boundaries. Therefore, before stacking 
the Pari data, such jumps were removed by subtract- 
ing the background from each quadrant using our home- 
grown routines. The quadrant offsets were reduced to 
~ le~ level after this process. 

The final mosaics were created following the standard 
drizzling procedures implemented in Multidrizzle, and 
were normalized to unit exposure time (count per sec- 
ond). A square kernel was used in drizzling, and the 
linear size of drop was set to 0.9 (the "finaLpixfrac" pa- 
rameter of Multidrizzle) . To get a finer spatial resolution, 
the drizzle scale (the "finaLscale" parameter) was set to 
0.6, which resulted in an output pixel scale of 0.03"/pixel. 
For each field, the output mosaics in different passbands 
were all registered to the same reference position in the 
process. The absolute astrometry of the final mosaics 
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was calibrated by using the compact objects visible in 
the Digital Sky Survey images. 

In order to provide a reference frame for Pari NIC3 
data reduction, another version of Pari ACS mosaics was 
also created by setting drizzle scale to 1.0, i.e., preserv- 
ing the native pixel scale (0.05'7pixel) of the ACS/WFC. 
We will refer to these lower-resolution mosaics and the 
higher-resolution ones mentioned above as "LOWRS" 
and "HIGHRS" ACS mosaics, respectively. They each 
serve different purposes in the analysis process. 

As the original pixels were resampled during drizzling 
(and during geometric distortion correction as well), the 
final mosaics all have correlated pixel noise. The signal- 
to-noise ratio (S/N) of a given object derived directly 
from the final mosaic would therefore be overestimated, 
and hence its photometric error would be underesti- 
mated. We used a routine, kindly provided to one of us 
(HY) by Dr. Mark Dickinson, to calculate the correla- 
tion amplitude and hence to calculate proper statistical 
error associated with each pixel. The error map (so- 
called "RMS" map) derived in this way are widely used, 
e.g., in the data analysis of GOODS (e.g., Dickinson et 
al. 2004). 

3.2. NIC3 data reduction 

All data obtained by NICMOS show a number of per- 
sistent anomalies that cannot currently be handled by the 
OTFR pipeUne. Therefore, the OTFR-processed NIC3 
data of Subl and Sub2 fetched from the archive were 
further reduced before stacking. 

An additional count-rate and wavelength depen- 
dent non-linearity of NICMOS was recently identi- 
fied by the STScI NICMOS team (NICMOS Instru- 
ment Science Reports 2006-001, 002 and 003; see also 
http://www.stsci.edu/hst/nicmos/performance /anoma- 
lies/nonlinearity.html). This non-linearity was removed 
from the data by using the routine provided by the team. 

The quadrant bias of NICMOS also has a stochastic 
behavior similar to what described above for ACS. As a 
result, a bias-corrected NICMOS image using the stan- 
dard bias reference file usually show "pedestals" from 
quadrant to quadrant. Such pedestals were removed by 
using the PESKY task included in the STSDAS package. 

All OTFR-processed NIC3 data have a persistent, ad- 
ditive pattern across entire field. This pattern is different 
in different passbands, but is rather stable with respect 
to time. To remove this pattern, a "pattern image" in 
each band was first created by stacking all the pedestal- 
corrected images in this band. The pattern image was 
then subtracted from each individual image. 

After all the above steps, the individual images were 
combined using the series of drizzle tasks included in the 
DITHER package of STSDAS in IRAF. Each input sci- 
ence image was weighted by a weight image calculated 
from its associated error array stored in the [ERR,2] 
extension of the OTFR-processed image. For the pur- 
pose of matched-aperture photometry that will be dis- 
cussed below, the LOWRS Pari ACS zgso-band mosaic 
was binned in 2 x 2 and then used as the reference. The 
NIC3 images were all mapped (i.e., rotated and regis- 
tered) to this binned ACS image during the drizzling 
process. The drizzle scale was set to 0.5, which resulted 
in a final resolution of 0.10"/pixel, i.e., the same as the 
binned LOWRS ACS image. The final mosaics were nor- 



malized to unit exposure time as well. 

The associated RMS maps of all these NIC3 mosaics 
were derived in the way similar to what described above 
for the ACS mosaics. 

3.3. ACS source detection 

As compared to the ACS mosaics used in YWC03, 

the new mosaics created through the imporved reduc- 
tion process (see §3.1) are of much higher quality, and 
hence we can push to a significantly lower thresold for 
source detection. 

SExtractor (Bertin & Arnouts 1996) was run on the 
HIGHRS mosaics in dual-image mode to extract source 
catalogs for both Pari and Par2. "MAG_AUTO" mag- 
nitudes were used to best approach the total magni- 
tudes. A number of source catalogs were obtained for 
different purposes. The zyys-band and ^gso-band mo- 
saics were used alternatively as the detection images to 
perform "matched- aperture photometry", i.e., measur- 
ing the magnitudes of a given source in different bands 
through the same aperture defined by its appearance 
on the detection image. We will refer to these catalogs 
as "z775-based" and "zsso-based" ACS catalogs, respec- 
tively. Source detection was performed at the threshold 
of 0.8 (J, using a 5 X 5 Gaussian filter with a full-width- 
at-half-maximum (FWHM) of 2.5 pixels. A minimum 
number of 4 connected pixels were required for a source 
to be included. 

As Pari was observed in parallel mode during the very 
early stage of ACS, the data were all taken at a high gain 
setting of 4.0 e^/ADU, probably for safety reason (i.e., 
avoiding saturation). In order to make the zcropoints 
derived at a gain setting of 1.0 e~/ADU universally ap- 
plicable, the ACS OTFR pipeline multiplies the flat-field 
reference file by the gain of the object image before flat- 
field correction. In the end, an OTFR-processed ACS 
image always has an effective gain of 1.0 e~/ADU. When 
YWC03 was written, unfortunately, the authors were not 
aware of this — at that time poorly-documented — step 
in the OTFR processing, and made a redundant correc- 
tion for the gain by adding 2.5 x loglQ{A.O) — 1.505 to 
the zeropoints. As a result, while the candidate selection 
in YWC03 remains valid, all magnitudes reported in that 
paper were too faint by 1.505 mag. This error has been 
reported in a footnote in YW04. 

In this current study, the zeropoints published in the 
latest ACS handbook were used; namely, 26.068, 25.654 
and 24.862 mag for 5475, 1775 and 2:850, respectively. The 
3 a limits within a circular aperture of 0.1 in radius are 
29.1 and 28.3 mag in ^775 and 2359 in Pari, and 28.9, 28.7 
and 27.9 mag in 5475, 4775 and 2:850 in Par2, respectively. 

3.4. NIC3 source detection 

Photometry of Subl and Sub2 NIC3 mosaics in Pari 
was also done in dual-image mode, using Jno and i?i6o 

mosaics alternatively as the detection images. Similarly, 
wc will refer to these catalogs as " Jno-bascd" and "iJif,o- 
bascd" catalogs, respectively. The detection filter was a 
5x5 Gaussian filter with a FWHM of 3 pixels, the de- 
tection threshold was set to 0.4 tr, and a minimum of 
2 connecting pixels were required. Again, MAG_AUTO 
magnitudes were adopted. As the magnitude zeropoints 
provided by the OTFR pipeline were no longer valid after 
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the non-linearity correction described in §3.2, wc followed 
the recipe given by the NICMOS team to derive the ze- 
ropoints as 23.242 and 23.139 mag in Jno and Hiqq, 
respectively. 

During this process, ACS magnitudes of the NIC3 
sources were also extracted in order to properly measure 
their ACS-to-NICMOS colors for later analysis. As the 
NIC3 mosaics were created by registering to the reference 
frame defined by the 2x2 binned LOWRS ACS mosaics, 
matched-aperture photometry between NIC3 and ACS 
passbands was fairly straightforward. The detection im- 
ages were still the NIC3 mosaics, and the object images 
were 2x2 binned ACS mosaics that were convolved by 
the PSF of the Jno images (derived using the TinyTim 
tool). The binning operation was sum, which preserves 
the original ACS magnitude zeropoints. 

4. OVERDENSITY OF 7-DROPOUTS 

Using the newly produced ACS mosaics of Pari, new 
NIC3 data obtained at two pointings within Pari, and 
a new ACS parallel field Par2 for comparison, here we 
will re-examine the overdensity of i-dropouts in Pari re- 
ported by YWC03. 

4.1. Overdensity of i- dropouts in Pari 

YWC03 reported 30 i-dropouts in Pari, whose zs5o 
magnitiidcs range from 26.8 to 28.3 mag. Using the 
correct zeropoints, their true zgso-band magnitude range 
should be 25.3 to 26.8 mag (see §3.3). The corresponding 
cumulative surface density thus is ~ 2.7 per arcmin^ to 
Z850 = 26.8 mag. Taken at its face value, this number is 
3 4 times higher than those in other deep, high galactic- 
latitude ACS fields, including the Hubble Ultra Deep 
Field (HUDF; Beckwith ct al. 2006). The i-dropouts 
in YWC03 were selected based on their invisibility in the 
2775 image, i.e., by selecting sources that have S/N > 5 
in Z850 but have S/N < 2 in 4775. As a result, about 2/3 
of those objects have 1775 — Zgso > 2.0 mag. Other stud- 
ies utilizing ACS data (including YW04) all have adopted 
less stringent selection criterion of 2775 — 2:850 > 1.3 mag 
or 2775 — 2850 > 1-5 mag. Therefore, the high surface 
density of YWC03 is even more difficult to be reconciled 
with all the other later work. 

To study this problem in detail, here we reselect i- 
dropouts in Pari from the new 2850-based ACS cat- 
alog as described in §3.3. We only consider sources 
with S/N > 5 in zsso, and adopt the color criterion of 
*775 ^ zs5o > 1.3 mag for case of comparison with other 
results. The candidates thus selected are visually exam- 
ined to eliminate false detections such as residual cosmic 
rays and image defects. The final t-dropout sample con- 
sists of 356 objects to 2350 = 28.0 mag. This is much 
deeper than what we achived in YWC03, and is the re- 
sult of the much improved data quality. In particular, 
198 of them have 2850 < 26.8 mag (among which 21 have 
^775 — ^850 > 2.0 mag), which is a factor of 20 more than 
the i-dropouts selected in the HUDF (YW04; Bunker et 
al. 2004) to the same depth. This confirms and re- 
inforces the unusually large overdensity of z-dropouts 
seen by YWC03. The distribution of these i-dropouts in 
Pari are shown in Fig. 2. 

A necessary condition for an i-dropout being a legiti- 
mate candidate of galaxy at 2 « 6 is that it must be a 
dropout in bluer bands. However, Pari does not have 



any data to the blue of 1775. Fortunately, as a subset 
of these i-dropouts have NIC3 observations, we can use 
their ACS-to-NIC3 colors to judge whether they are con- 
sistent with being galaxies at 2 w 6. 

4.2. Optical-to-IR colors of i-dropouts in Pari 

In total, 21 i-dropouts are detected in our NIC3 data. 
Fig. 3 displays the image cut-outs of one such ob- 
jects. The optical-to-IR color-color diagrams of all these 
sources are shown in the left panel of Fig. 4. For com- 
parison, the right panel of Fig. 4 shows the similar color- 
color diagrams from the HUDF. 

It is clear that the i-dropouts in Pari reside in the 
similar regions in the optical-to-IR color space as the i- 
dropouts in the HUDF. In particular, they are well sepa- 
rated from brown dwarf stars and E/SO galaxies at lower 
redshifts, which are the most common contaminators to 
2 w 6 galaxy samples selected by the drop-out technique. 

4.3. A similar overdensity of i-dropouts in Par2 

In order to compare with the results obtained in Pari, 
we perform z-dropout selection in Par2 in a similar way. 

As this field has 9475 data, an additional constraint of 
non-detection in 17475 at 2 cr level is imposed for the se- 
lection. The final i-dropout sample consists of 20 objects 
to 2850 = 26.5 mag, which, while much less than the ones 
in Pari, is still a factor of ~ 7 more than those found in 
the HUDF to the same depth; in other words, the over- 
density of i-dropouts also seems to present in Par2. 

5. ALTERNATIVE INTERPRETATION OF THE 
OVERDENSITY 

The analysis in §4 indicates that the i-dropouts in Pari 
are consistent with being galaxies at z rs 6, and that the 
overdensity is also present in Par2. If this interpreta- 
tion is true, we will have to draw the conclusion that 
the overdensity of i-dropouts in this region is due to an 
extremely unusual concentration of 2 « 6 galaxies over 
~ 10 Mpc scale. In this section, however, we present 
another possibility. 

5.1. Excess of field population in Pari & Par2 

The alternative interpretation is triggered by an un- 
usual field population seen in both Pari and Par 2. This is 
illustrated in Fig. 5, which shows the 2850 vs. i775— 2850 
color-magnitude diagrams in these fields using the 2859- 
based catalogs, together with the corresponding dia- 
grams constructed in the HUDF based on the photome- 
try from YW04. The left panel only plots the sources at 
S/N > 5, while the right panel includes all the sources 
at S/N > 3. A striking feature in this figure is the large 
concentration of objects at the faint end in the two par- 
allel fields, which is not seen in the HUDF (and other 
deep ACS fields). Due to the shallower depth in Par2, 
this population is at the verge of being detectable at 
S/N > 5; however, it gets very prominent if the thresh- 
old is lowered to S/N > 3. YWC03 noticed only a hint 
of this excess in Pari; the soiirce detection back tlicni did 
not reach a sufficient depth that could unambiguously 
confirm its existence. 

This excess of faint field objects in Pari and Par2 is 
present in all passbands. To quantify this excess. Fig. 
6 compares the surface densities of the detected sources 
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in these two fields and the ones in the HUDF, derived 
at two different threshold levels oi S/N > 3 (left panel) 
and S/N > 5 (right panel). Again, while the excess in 
Par2 is only detectable when the sources of low S/N are 
included, the excess in Pari is obvious even at a high 
source detection threshold. This excess starts at about 
26.0 mag in both Pari and Par2, and extends to a much 
fainter level. In particular, the deep mosaics in Pari 
allow this excess to be traced to at least ~ 28.2 and 27.2 
mag in ^775 and 2:850, respectively. To further ensure 
the existence of the excess, wc split the available images 
of Pari into two groups and created shallower stacks in 
both bands. The same excess was detected as well. 

The i-dropout surface densities are also superposed in 
Fig. 6. We point out that the i-dropout excess in the two 
parallel fields resides in the region where the field popu- 
lation excess happens. This is further discussed later in 
the paper. 

5.2. Resolving diffuse light from interacting galaxy pair 

What is the nature of this field population excess? The 
answer probably lies in the proximity of Pari and Par2 

to the galaxy pair M60/NGC4647, which is several ar- 
cminutes away to the South (sec Fig. 1). Extended dif- 
fuse emissions of different natures have been detected 
around some nearby galaxies, for example in NGC 5907 
(Shang et al. 1998; Zheng et al. 1999) and in M31 (e.g., 
Ferguson et al. 2002). While there is no previous re- 
port in the literature, we speculate that such extended 
diffuse emission also exists around M60/NGC4647, and 
propose that the excess population is due to the fact 
that we have resolved the diffuse light extending from 
this interacting pair. If we adopt an average distance 
of ~ 15.7 Mpc to M60/NGC4647 (from NASA/IPAC 
Extragalactic Database). Pari and Par2 arc ^ 32 and 
69 kpc away from the center of pair, respectively. There- 
fore, what we have detected are most likely the halo stars 
of the two galaxies and/or the tidal "debris" repelled 
from this pair during the course of the interaction. 

Here we exam the property of this excess more closely. 
We only use Pari for this purpose, as the data in this field 
are much deeper than in Par2. Fig. 7 plots the surface 
densities of S/N > 5 sources in Pari as functions of 
magnitudes in 2775 (left) and zgso (right), and compares 
them with the results in the HUDF. The insets show the 
net excess after subtracting the surface density calculated 
in the HUDF. Tab. 1 lists the surface densities of the net 
excess as a function of magnitudes. 

Taken at face value, the net excess in Pari as shown in 
Fig. 7 has a total magnitude of 15.0 and 15.7 mag in 1775 
and 2850, respectively. The average surface brightness 
from the net excess (i.e., spreading the total integrated 
light of the net excess over the entire field of Pari), on 
the other hand, is 26.6 and 27.2 mag/arcsec^ in 2775, 
and 2850, respectively. If this excess is isotropic with re- 
spect to M60/NGC4647, the total magnitude integrated 
over an annulus 3.4 in width (i.e., about the width of 
one ACS field) at the distance of Pari (centering on the 
galaxy pair) is 12.2 mag in 2775 and 12.9 mag in 2850- 
For comparison, the total flux of M60 in Johnson I-band 
is 8.46 mag (or 0.925 Jy, i.e., 8.98 mag in AB) within 
84 aperture (Boroson, Strom & Strom 1983). There- 
fore, the total diffuse light within this annulus is about 



a few percent of the central galaxies. 

The net excess peaks at about 27.1 and 26.6 mag in 4775 
and 2850. If these objects are indeed at the distance of 
Virgo Cluster, these values correspond to absolute mag- 
nitudes of —3.9 and —4.4 mag, respectively, which are in 
the regime of giants/supergiants or low luminosity glob- 
ular clusters. 

We shall point out that the excess that we detected re- 
lates more closely to M60/NGC4546 rather than to the 
Virgo Cluster in general. A supporting evidence comes 
from the spatial gradient of this excess, which decreases 
in amplitude when moving away from the pair. This is 
demonstrated in Fig. 8, which plots the spatial distribu- 
tion of all S/N > 5 sources at 2350 > 26.0 mag detected 
in the 2850-based catalog of Pari. The gradient is obvi- 
ous: there are more sources to the south, which is the 
direction where the galaxy pair lies. 

Another piece of evidence is that the amplitude of this 
excess (see Fig. 7 and Tab. 1) cannot be attributed 
to the intracluster stellar populations of the Virgo Clus- 
ter alone. We only consider Virgo's intracluster glob- 
ular clusters (IGC) and red giants (IRG), as they are 
the only two type of stellar objects that arc within the 
luminosity range of the detected excess. The total contri- 
bution from Virgo's intracluster globular clusters (IGC) 
at > 26.0 mag, calculated using the Virgo's IGC surface 
density at the bright-end (e.g., Williams ct al. 2007a; 
Cohen et al. 2003) and the symmetry of globular cluster 
luminosity function, is only ^ 0.3 3.0 arcmin"^ (for any 
reasonable color transformation from F814W to F775W), 
depending on the proximity to a member galaxy. On the 
other hand, the deep HST surveys of Virgo's IRG (e.g., 
Ferguson, Tanvir & von Hippel 1998; Durrell et al. 2002; 
Williams et al. 2007b) show that their cumulative surface 
density to ~ 28.0 mag can be as high as ~ 110 arcmin"^ 
(again for any reasonable F814W to F775W transforma- 
tion). However, this is still a factor of 5-6 lower than 
what shown in Fig. 7 and Tab. 1. 

5.3. Relation between i-dropout overdensity and field 
object overdensity 

We suggest that the i-dropout overdensity observed in 
Pari & Par2 might be related to the field object overden- 
sity in the same regions. The field object excess becomes 
significant at 26.0 mag in both bands, and extends to 
at least 28.2 mag in 2775 and 27.2 mag in 2850 in the 
deeper Pari images (sec the insets in Fig. 8) where the 
excess can be traced to a faint level at high confidence. 
The i-dropouts selected in these fields resides in the 2850 
magnitude range that the field object excess spans (Fig. 
6). The spatial distribution of the i-dropouts also seems 
to follow the gradient that seen in the distribution of the 
excess field population (see Fig. 2 & 7). 

This excess field population is also red in color. Fig. 
9 shows the surface densities of the sources in the above 
mentioned magnitude ranges as a function of ^775— 2850 
color. For comparison, the surface densities of the HUDF 
sources in the same magnitude ranges are also shown. 
The insets plot the net excess using the HUDF re- 
sult as the reference. The excess population peaks at 
^775— 2:85o~ 0.35 mag in the 1775 image, and peaks at 
*775— Z850~ 0.60 mag in the 2350 image. Note that the 
distribution extends continuously beyond 1775— 2350 > 1.3 
mag, which is the regime where the i-dropouts are se- 
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lected. 

If the i-dropouts are indeed some stellar objects re- 
lated to M60/NGC4647, what type of stars could they 
be? Thus far we have not yet been able to identify the 
satisfactory candidates. As Fig. 4 shows, brown dwarfs 
are not likely to contribute much to our i-dropout sam- 
ple because their optical-to-IR colors are different. In 
addition, their luminosities are too low. The only pos- 
sibility then left is dusty stellar population, for exam- 
ple, Asymptotic Giant Branch (AGB) stars. However, 
the optical-to-IR colors of AGB stars are quite different 
from the i-dropouts as well. This is illustrated in Fig. 
10, where the synthesized colors of AGB stars generated 
from the model spectra of Lancon & Mouhcine (2002) are 
superposed on the color-color diagrams of the i-dropouts 
in Pari. We should also point out that no external dust 
extinction to any known types of stars could produce 
such colors; no know extinction law could create such 
a red color in 1775 — ^§50 and yet keep zs5o — Jno and 
Jiio — HiQo colors around zero. 

The red, excess field objects that we have discovered 
here seem to have some connection to the so-called "Red 
Halo Phenomenon" (Zackrisson et al. 2006), which re- 
ferrs to the very red colors of the halos found around 
some galaxies (e.g., Zibetti, White & Brinkmann 2004; 
Zibetti & Ferguson 2004; Taylor et al. 2005). Such red 
halos cannot be easily explained by standard stellar pop- 
ulations. The Z775 — Z850 color shown in Fig. 9 quan- 
titively agrees with that of the halo around an edge-on 
disk galaxy in the HUDF reported by Zibetti & Fergu- 
son (2004). Therefore, the excess field objects that we 
see might represent the first example where such red ha- 
los have been resolved, and the i-dropouts might be the 
most extreme sources among this population that is most 
difiicult to explain. 

6. SUMMARY 

In this paper, we use new and existing HST data to 
closely examine the large overdensity of i-dropouts that 
YWC03 found in a deep ACS parallel field. The new 
data include NIC3 imaging of a subset of the i-dropouts, 
and a set of new ACS parallel data in a nearby region 
obtained during the NIC3 observations. We confirm that 
a factor of ^ 20 overdensity of i-dropouts does exist in 
the parallel field of YWC03, and that the optical-to-IR 
colors of these objects are similar to those selected in 
other deep ACS fields such as the HUDF. 

However, these z-dropouts might not be galaxies at 
2; ~ 6 as one would expect. During the course of our 



investigation, we have found a large excess of discrete 
sources in both the field of YWC03 and the new ACS 
parallel field. This excess is most likely due to a stellar 
population that is related to M60/NGC4647, an interact- 
ing galaxy pair in the Virgo Cluster, which lies several 
arcminutes to the south of the ACS parallel fields. At the 
distance of the Virgo Cluster, the ACS parallel fields arc 
several tens of kpc away from the galaxy pair, and hence 
the stellar objects are likely tidal "debris" expelled from 
the pair and/or distant halo stars belong to these two 
galaxies. Both the spatial distribution and the ampli- 
tude of the excess indicate that it is a population related 
more closely to the galaxy pair rather than to the Virgo's 
intracluster populations. 

The i-dropouts in the field of YWC03 are within the 
magnitude range where the field object excess occurs. 
Their distribution also shows similar spatial gradient of 
the field object excess. Most importantly, the field object 
excess is red in color, and the red wing of its color distri- 
bution extends well into the regime where the i-dropouts 
were selected. For these reasons, we believe these i- 
dropouts are part of the stellar objects that constitute 
the field object excess rather than galaxies at z ~ 6. 

Nevertheless, we should caution that we still cannot 
completely rule out the possibility that these i-dropouts 
might indeed he z ^ 6 galaxies. On one hand, we do 
not have any definite evidence against this interpretation. 
On the other hand, we have not yet found a stellar pop- 
ulation that can reasonably explain the colors of these i- 
dropouts. To unambiguously identify their nature, spec- 
troscopic observations of these objects are in demand. 
As the brightest sources in our sample arc brighter than 
26.0 mag in zs5o, this might be feasible with the existing 
spectrographs at our current 8-lOm class telescopes. 
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TABLE 1 

Surface density of the S/N > 5 source net excess in Par! as a function of magnitude 



magnitude 


S(i775) (arcmin 


^{zs5o) (arcmin 2) 


25.55 


0.51 


1.35 


25.65 


1.96 


0.93 


25.75 


2.52 


2.62 


25.85 


5.13 


2.51 


25.95 


5.54 


4.67 


26.05 


7.11 


6.99 


26.15 


9.60 


13.13 


26.25 


14.58 


15.86 


26.35 


18.93 


19.76 


26.45 


29.74 


29.16 


26.55 


34.35 


30.55 


26.65 


36.93 


24.61 


26.75 


38.54 


18.89 


26.85 


42.51 


12.57 


26.95 


43.04 


8.78 


27.05 


47.17 


7.04 


27.15 


45.84 


8.32 


27.25 


43.59 


2.24 


27.35 


36.64 




27.45 


39.49 




27.55 


37.16 




27.65 


26.10 




27.75 


24.45 




27.85 


18.44 




27.95 


11.06 




28.05 


9.66 




28.15 


0.58 





■ The net excess is as shown in the insets in Fig. 7, which is calculated using the count in the HUDF as the reference. 
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Fig. 1. — The locations of the fields studied in this work are shown on top of the image from the Digital Sky Survey. Both fields are 
several arcminutes away from the interacting galaxy pair, M60 & NGC4647, which are members of the Virgo Chister. Pari was taken in 
parallel mode by the ACS when the WFPC2 was observing NGC4647, while Par2 was taken by the ACS in parallel when NICMOS was 
observing two regions in Pari. 
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Fig. 2. — The left panel shows the locations of the j-dropouts in Pari as red dots on top of the zgso-band mosaic of this field. The yellow 
dots in the right panel show the locations of the i-dropouts that are detected the NIC3 mosaic, which is the image in the background (note 
that the objects that are seemingly close to the edge are actually well within the region of good coverage). 
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Fig. 4. — The optical-to-IR colors of the i-dropouts in Pari (left) are compared with those of the i-dropouts in the HUDF (right). The 
i-dropouts are the red points with error bars, while the field objects are shown as the black dots. M-, L- and T-type brown dwarf stars 
are the green, blue and brown asterisks, respectively. The colors of E/SO galaxy at z = 1-3 are shown as the yellow crosses. Two model 
galaxies at 5.5 < z < 6.5, constructed using the models of Bruzual & Chariot (2003) of different star formation histories, are shown as the 
open magenta squares. The j-dropouts in Pari occupy a similar region in the color space as do those in the HUDF. 
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Fig. 5. — Comparison of the zsso vs. 2775—2:850 diagrams between the parallel fields and the HUDF. The objcets plotted in the left 
panel are limited to the sourees with S/N > 5, while those in the right panel are limited to the sources with S/N > 3. The dotted lines 
show the •i775— -2850= 1-3 eolor threshold above whieh the i-dropouts are selected. A unique feature in the two parallel fields is the large 
excess objects at the faint end, which is not seen in the HUDF. Due to its shallower depth, Par2 does not show this feature very clearly if 
only the sources at S/N > 5 are included; the excess becomes prominent when sources to S/N = 3 are also included. 
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Fig. 6. — Comparison of the field source surface densities in the two parallel fields (solid and dashed black curves) and in the HUDF 
(blue dotted curves) at two threshold levels. The irrs-based catalogs are used for the statistics in iyrs-band (top), while the zgso-based 
catalogs are used for the statistics in zgso-band (bottom). The left panel shows the statistics based on sources with S/N > 3, while the 
right panel shows the statistics based on sources with S/N > 5. The object excess in the two parallel fields is obvious. The surface densities 
of i-dropouts (all have S/N > 5 in 2350) in Pari, Par2 and HUDF are also shown as the magenta, cyan and red curves, respectively. 
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Fig. 7. — Field source excess of Pari as a function of magnitude. The solid histograms are the surface densities in Pari, and the dotted 
histograms are those in the HUDF. The i775-bascd catalogs arc used for the figure in the left panel, while the zsso-based catalogs are used 
for the one in the right. Only sources at S/N > 5 are included for the statistics. The insets show the net excess using the HUDF as the 
reference. The excess extends to ~ 28.2 and 27.2 mag in in^ and zgsOi respectively. The average surface brightness of the net excess 
corresponds to 26.6 and 27.2 mag/arcsec'^ in 1775, and zgso respectively. The total magnitudes of the net excess integrated over Pari are 
1775=15.0 and 2850=15.7 mag. 
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Fig. 8. — The spatial location of the excess field population has a clear gradient towards the direction of the interacting galaxy pair 
M60/NGC4647. This figure plots the positions of the Zs5()-based sources with S/N > 5 and Zfi-,Q> 26.0 inag, which is the brightness level 
where the excess starts to be significant. North is up and East is to left. M60/NGC4647 is to the south of the field (see also Fig. 1). Note 
that the gap in the middle is caused by the gap between the two CCD chips of the ACS/WFC. 
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Fig. 9. — The excess population in the parallel fields is red in color. The solid histograms show the iTTs— Z850 color distributions of Pari 
sources that are within the magnitude ranges where the excess occurs, while the dotted histograms show the distribution of the HUDF 
sources within the same magnitude ranges. The insets show the distributions of the net excess using the HUDF as the reference. The 
*775-based catalogs are used in the left panel, and the zggo-based catalogs are used in the right panel. The red wing of the excess population 
in 2850 continuously extends to the regime where the i-dropouts are selected. 
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Fig. 10. — While the i-dropouts in the two parallel fields might well be part of the excess stellar population related to M60/NGC4647 
instead of galaxies at z ^ 6, we have not yet been able to identify the type of stars that can create the colors of the i-dropouts. The most 
promising candidates are AGB stars, which have the necessary luminosities and are red in colors. However, AGB stars still cannot fully 
explain the colors of the i-dropouts. To illustrate this, the same color-color diagrams as in the left panel of Fig. 4 are shown here again, 
and the colors of AGB stars (cyan squares) from the models of Lancon & Mouchine (2002) arc superposed. To avoid confusion, only the 
models that satisfy irrs— 2850> 1-3 mag are superposed in the right panel. 



